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Abstract 

We investigate the correlations of initially separable probability distributions in a globally pure 
bipartite system with two degrees of freedom for classical and quantum systems. A classical version 
of the quantum linear mutual information is introduced and the two quantities are compared for a 
system of oscillators coupled with both linear and non-linear interactions. The classical correlations 
help to understand how much of the quantum loss of purity are due to intrinsic quantum effects and 
how much is related to the probabilistic character of the initial states, a characteristic shared by 
both the classical and quantum pictures. Our examples show that, for initially localized Gaussian 
states, the classical statistical mutual linear entropy follows its quantum counterpart for short 
times. For non-Gaussian states the behavior of the classical and quantum measures of information 
are still qualitatively similar, although the fingerprints of the non-classical nature of the initial 
state can be observed in their different amplitudes of oscillation. 
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I. INTRODUCTION 



EntaiiErlement has been a focus of intense investigation in the recent years due to its 
relevance in quantum computation and quantum information From a 

fundamental point of view, entanglement is considered "the characteristic trait of quan- 
tum mechanics, the one that enforces its entire departure from classical lines of thought" 
(Schrodinger) 7]. Entanglement reflects quantum correlations in the Hilbert space, where 
state vectors are associated with probability distributions. In classical mechanics, on the 
other hand, states are described by points in phase space, and no intrinsic probabilistic 
character exists. 

Probability distributions can be introduced in classical mechanics with the Liouville for- 
malism, and statistical averages on ensembles can be calculated. Initially independent prob- 
ability distributions can become correlated when evolved by a classical Hamiltonian with 
interaction terms. These correlations represent the dynamical emergence of conditional 
probabilities in the classical statistical world. 

In this paper we consider the time evolution of a bipartite system whose Hilbert space 
is the direct product of two subspaces. Initial states that are the tensor product of kets in 
each subspace, usually evolve to non-product states, generating coherences and entangling 
the two subsystems. A measure of the non-separability between them is provided either by 
the Von Neumann or the linear entropies. In the case of a globally pure bipartite system 
we argue that the latter is a convenient quantum measure of non-separability. Next, as the 
main point of our work, we consider the time evolution of probability densities in the phase 
space of the corresponding bipartite classical analogue. Inspired on the quantum problem 
we propose a measure of the classical correlations generated by the statistical ensembles of 
classical trajectories. We want to measure the non-separability of the classical probability 
distribution evolving in time via Liouville equations. Moreover we want to compare behavior 
of this classical measure with the loss of purity of the quantum system, particularly at short 
times. This comparison helps to identify the time at which intrinsic quantum effects, such 
as interferences, begin to be important. 

In our approach we first define a classical statistical quantity and compare its behavior 
with the quantum linear entropy (QLE). Due to its formal similarity with the quantum linear 
entropy, we call it classical statistical linear entropy (CSLE). A similar classical quantity has 
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recently been used in the study of quantum-classical correspondence of intrinsic decoherence 
^. In spite of the direct correspondence between the QLE and the CSLE, the latter may 
not be symmetric between the two subsystems that make up the global bipartite system, 
i.e., the CSLE of each subspace does not contain in general the same information. Also, the 
CSLE can assume negative values and is restricted to initially separable distributions, not 
being a good measure in the case of distributions that start off non-separable. We therefore 
define a more general classical measure, the classical statistical linear mutual information 
(CSLMI), based on the concept of quantum mutual information (QMI) [9j, that is always 
symmetric and positive, and thus, best suited to measure the classical separability. We show 
that for two bi-linearly coupled harmonic oscillators with Gaussian initial distributions, the 
classical and quantum results coincide. For non-linear coupling the classical and quantum 
entropies remain close for short times for both regular and chaotic regimes. 

The outline of this work is as follows. In Section II. A we give some general definitions and 
show that the quantum linear entropy can be expressed in terms of the non-diagonal terms of 
the density operator, providing a good measure of the entanglement between the subsystems. 
In Section II. B we present our definitions of the classical statistical linear entropy with some 
considerations about the initial probability distributions. Next we introduce the quantum 
linear mutual information and its classical statistical analogue. This allows us to treat cases 
where the two subsystems have different types of initial distributions in phase space. Section 
III is reserved to the comparison between the QLMI and the CSLMI for a system of two 
oscillators coupled with the following types of interactions: (1) bi-linear; (2) non-linear ; (3) 
bilinear within the rotating wave approximation (RWA). As initial states we consider both 
Gaussian and non-Gaussian distributions. Finally in Section IV we present some conclusions 
and discuss the adequacy and limitations of the present approach. 

II. THEORY 

A. Quantum Linear Entropy 

The density operator corresponding to an initially (normalized) pure state l^po) is 

p{t) = u{t)\^Po){Hu\t), (1) 
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where U{t) = e~'^*/'* is the evolution operator and H is the Hamiltonian. For a globally 
pure state, there are constraints connecting the diagonal and off-diagonal matrix elements of 
the density operator. This follows from the property Trp^ = Trp = 1, which can be written 
explicitly as 

$ZP-' + $Z5ZlPifcl' = $ZP- = l' (2) 

i i k^i i 

where we have separated the diagonal and off-diagonal terms in the left hand side. From 
this relation we define 

i k^i i 

Consider now a bipartite system whose Hilbert space is the direct product of two sub- 
spaces: £ = £i® £2. A measure of non-separability between these subsystems for a globally 
pure state is provided either by the Von Neumann or by the linear entropies 0, Q|. The 
partial trace of p on system 2 defines the operator 

Pi{t)=tr2{p{t)) . (4) 

The corresponding QLE is given by 

S,{t) = l-tT,{pl{t)). (5) 

The operator P2{t) and the QLE 5*2 (t) are similarly defined. 

For if = ifi ® I2 + li ® i/2 both pi(t) and p2(t) are projectors and Si{t) = S2(t) = 0. If 
couplings between the sub-spaces exist, they force the global state to evolve to non-product 
states, entangling the subsystems and producing non-zero values for Si(t) and 5*2 (t). 

In fact, for a globally pure bipartite system Si = S2 = d{p). To see this we write 

\^) = J2^M^)\b^), (6) 

i 

where {|aj)} and are orthonormal basis of subsystems 1 and 2 respectively, and {Aj} 

is the Schmidt spectrum Then, 

P = ^PiM)h){(^3\{bj\, (7) 

id 
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with pij = AjA*. The reduced density of system 1 becomes 

Pi = ^Pii\ai){ai\, (8) 

i 

and similarly for system 2. The subsystem linear entropy becomes trivial in this basis: 

S, = l-Tn{pl) = l-J2pl (9) 

i 

with a similar result for 5*2. 

Comparing this result with Eq.Q we see that 

5l = 52 = 5^5^|p.fcP = f/(p). (10) 

Therefore, as long as any globally pure state of the system can be decomposed in the 
Schmidt basis, we can conclude that the linear entropy indeed measures coherences between 
the subsystems. Furthermore, because of the total conservation of coherences (j21), the linear 
entropy seems to be the most natural measure in the present situation. 



B. Classical Systems 

A measure of non-separability for classical systems can be defined only at a statistical 
level by considering ensembles of initial conditions Q, Following Wehrl jl5| we 
define a quantity that we call classical statistical linear entropy. 

Consider a system with two degrees of freedom described by the classical Hamiltonian 
function Ti. Consider also several copies of this system with initial conditions distributed 
according to the ensemble probability distribution P{x,t = 0), where x = {qi,pi,q2,p2)- 
The classical time evolution of P{x, 0) is obtained via Liouville's equation 

f^{«.n (11) 

whose solution is 

P{x,t) = P{<P;\x),0), (12) 

where 0jr^(a;) = Xq is the initial condition that propagates to x in the time t and (f>t is the 
phase space flux, so that x = (f)t{xo)- In words, the numerical value of the probability P at 
point X and time t has the same numerical value of the probability at point Xq of the initial 



distribution. This value is carried over to x in the time t. Liouville's theorem guarantees the 
conservation of probabihty at all times. For integrable systems (f)t{xo) might be determined 
analytically, otherwise numerical calculations can be performed. Notice that P{x, t) is itself 
a constant of motion, but this is not enough to guarantee the integrability of the system, 
since it is not generally in involution with the Hamiltonian. 
The marginal probability distribution 

Piiqi,Pi,t) = J dq2dp2P{qi,Pi,q2,P2,t) (13) 

allows us to define the classical statistical linear entropy 

qd{+\ - 1 _ j dqidpiP^{qi,pi,t) 

! dq,dp,P^q,,p,,Q) ■ ^'^> 

The normalization is necessary for dimensional reasons and to guarantee that 5*^(0) = 0. 

The initial classical distribution corresponding to a given quantum pure state p(0) is 

chosen as the coherent states phase space projection 

P{x, 0) = a2|p(0)|a2, Oil), (15) 



where iV is a normalization constant and = (g^ + tpi)/\/2h is the usual complex 
parametrization of the coherent states. Eq. (fT3j) is the normalized Husimi distribution, 
which is positive definite by construction. It is known that the Husimi distribution does 
not reproduce the correct marginal probabilities. However, the constraint of a positive 
probability distribution excludes the Wigner function as a classical initial distribution. 

Eqs. (jn} and show that the quantum linear entropies 5*1 and 5*2 can be obtained only 
from the diagonal elements of the global density matrix. Diagonal elements, on the other 
hand, have classical analogues. Therefore it makes sense to compare the dynamics of 5*1 
and 5*2, which may be written either in terms of only off-diagonal or only diagonal elements, 
with the classical correlations. 

An important difference between the classical and quantum linear entropies can be made 
explicit in the coherent state basis. For one dimensional systems we obtain 

S{t) = l- / —\{a\U{t)p\0)U^{t)\a)]. (16) 



On the other hand, defining the Liouvillian operator C such that {7i, P} = CP, Eq. ljllll can 
be formally integrated and the CSLE can be written as 



S'^{t) 



d'^a 

TT 



{e^'{a\pma)y 
M 



(17) 
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where M = N'^/ir. In particular, at t = the integrand of the classical entropy depends on 
{a\p{0)\a)'^ instead of {a\p'^{0)\a). Another important difference is of course in the dynamics: 
the quantum evolution is determined by non-commuting operators, which brings a number 
of corrections to the classical formalism. For more degrees of freedom, although the time 
evolution cannot be expressed in such a simple way, the two differences pointed out above 
remain true. 

C. Quantum and Classical Linear Mutual Information 

Although the CSLE seems to be the natural classical analague of the quantum linear 
entropy, it is not symmetric between the subsystems. Indeed, if the initial probability dis- 
tributions of each subsystem are not equal (for example a Gaussian distribution for one 
subsystem and a Poisonian distribution for the other) we find that Sf{t) ^ S2{t). The 
quantum linear entropy, on the other hand, always satisfies Si{t) = 5*2 (t). Another draw- 
back of the definition Eq. p4j) is that it gives 5*^(0) = 5*2' (0) = for any initial classical 
distribution, correlated or not. Thus, it is interesting to define a more general quantity that 
avoids these difficulties and that takes into account the contributions of both subsystems 
entropies symmetrically. At the quantum level we define the quantum linear mutual infor- 
mation (QLMI) that depends on the OLE as / = S{pi ® P2) — S{p). This is based on the 
Von Neumann mutual information j^Q, Ivn = Svn{pi ® P2) — Svn{p), where Svn is the 
Von Neumann entropy. By using Eq.(0) we can write the QLMI as (notice that the linear 
entropy is not additive ^IQj ) 

lit) = S,{t) + S2it) - Slit) S2it) - S{t), (18) 

where Si and S are the subsystem and the global linear entropies respectively. For pure 
initial states S{t) = and I{t) = Si + S2 — SiS2- 

We also define a quantity that we call classical statistical linear mutual information 
(CSLMI) as 

r\t) = Stit) + S^\t) - St\t) Sl^^it) - S'\t). (19) 

Once again S'^\t) = for pure states. Note that the above definition can also be written in 
form = (Pi P2) - S''\P) where ^"'(P) is given by Eq.^H} with Pi replaced by P. 
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We emphasize that the quantum and classical linear mutual information defined above 
are also non-separability measures, like the QLE and CSLE. This can be made explicit by 
re- writing Eqs.(fTS|) and (fTUI) as 

= 1- triplet) 0p2^it)l (20) 



r\t) 



JdxP^iO) 
J dxP,\t)P,\t) 

JdxP^iO) ■ ^ ^ 



The last equality is true since / dxP'^{0) = f dxP'^it). When the two subsystems have the 
same type of initial distributions, both I and P^ have the same contents of their respective 
linear entropies. Moreover, since P^ is symmetric with respect to the two subsystems, it 
does not present difficulties when the initial distributions are different. 
In the next section we compare I{t) with P\t) for three different cases. 



III. RESULTS 



In order to show that our definition of a CSLMI is physically sensible, we consider the 
following classical Hamiltonian with two degrees of freedom 

n{quPuq2,P2) = ni + n2 + ni , (22) 

where Hi = ^{Pi + ^ili) are harmonic oscillators and Hj is an interaction term. In what 
follows we shall calculate both CSLMI and QLMI for three different couplings: a simple 
bilinear, a non-integrable and a rotating- wave approximation (RWA). In the first and last 
cases the calculations can be performed analytically. 

1) Bilinear Coupling Hi{qi,pi,q2,P2) = M1Q.2 [ISf 
In this case, the classical equations of motion can be easily integrated and we find 

= giC+(t) + q^C.it) +p,S^{t)+p,S.{t), (23) 
Pi(t) = piC+(t) +p2C_(t) - q,S+{t) - q^S^t), 
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and similar expressions for q2(t) and p2it). C± and S± are quasi-periodic functions of time 

C±(t) = - (cos Qxt ± cos Qyt) , 
S±{t) = -^{Qx sin ± sin Qyt) , (24) 

where il^ = y/uj'^ + ^ and fiy = y/uj"^ — A. 

As initial phase space distribution we choose a Gaussian centered at Xc = {qic, Pic, i?2c, P2c), 

This is the classical density as given by Eq. (jl5p . corresponding to a coherent state initial 
wave-function. 

The probability distribution at time t is obtained using Eqs. (fT^ and (f^. The result 
can be written in the form 

P(x,t) = ^e-^"^^+2^^+^, (26) 

where the matrix A and the vector B are functions of C±, S± and h. Replacing Eq.()26|) 
into (fT^ we get 

P'{t) = 1 - ^ 2 ^ , (27) 

64/1^ det (/?)' det {a - 7/^-17^) ^ ^ 

where a, /3 and 7 are 2x2 matrices whose elements are combinations of C±, S± and T±, 
with 

/ sin , sin rij,t 

It is interesting to note that all the dependence on the center position of the initial distri- 
bution has gone. 

The reduced QLMI can also be analytically computed. Here, we simply write down the 
result: 



m = 1 



(28) 



The coefficients Di{t) are periodic functions of frequency fij. O is an 8 x 8 quasi-periodic 
matrix depending on both Vti and ^2- Expressions ()27p and ()28|) are actually identical, which 
shows that our definition of and the choice of its normalization are both appropriate. 
Fig. 1 displays an example for A = 0.9 that shows /i(t) = ifit)- 
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2) Nonlinear Coupling 

The exact coincidence between the classical statistical and quantum linear mutual informa- 
tions just presented, certainly has to do with both the quadratic nature of the Hamiltonian 
and the Gaussian initial distributions. To study the role of non-linearities we consider the 
interaction Hamiltonian 

'Hi ill, Pi, 12, P2) = -qiPiP2 + ^gi^g2^- (29) 

The total Hamiltonian 7i is a canonically transformed version of a well studied system known 
as the Nelson potential . Fig. El shows a typical mixed Poincare section for parameter 
values LUi = VOA, uj2 = and energy E = 0.05. 

The calculation of lf{t) (Eq.^) now has to be performed numerically. We use the 
same initial distribution, Eq. ^Kjl . xq = (pt^i^) was calculated using standard Runge-Kutta 
routines. The integrations in Eqs. pHj) and (fT^ was done both by Monte Carlo and by 
direct trapezoidal techniques; the results of the two methods agree in the time intervals we 
have considered. In Fig.(j2K) we show I{t) and the corresponding P\t). The initial density 
matrix is the direct product of two coherent states, and the initial classical probability 
distribution is that given by (fT3j) . The center of the coherent states is in the chaotic region 
of the corresponding Poincare section. The resemblance between the QLMI and CSLMI is 
quite good even after two oscillations. Fig(jnb) shows a similar calculation with the coherent 
states centered at the regular region. Once again the classical and quantum results coincide 
for short times. Surprisingly, the two entropy-like quantities agree for longer times in the 
chaotic case. For both the regular and chaotic cases, the classical and quantum calculations 
agree very well for short times, although the classical mutual information is systematically 
larger than its quantum counterpart for times larger than about 1. Also, the linear mutual 
information grows faster for the chaotic case than for the regular case, in accordance with 
similar previous results for the linear entropy Q]. Finally we note that the classica]_ linear 
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entropies Eq. lfT^ may become negative. Similar behaviors were reported in refs. 
The linear mutual information, on the other hand, is always positive and better suited for 
measuring the classical loss of separability. 

3) RWA Coupling: Hi{qi,pi,q2,P2) = A(gig2 + P1P2) 
This is the classical version of a rotating-wave approximation of the interaction Hamiltonian 
treated in the example (1). In this case, Gaussian distributions in each subspace evolve 
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coherently: Ii{t) = lf{t) = 0. However, this is a very particular case of a preferred basis 



state 



22| . and the same kind of coherent evolution is not expected for more general 



initial states. For instance, in the case of Fock states {\iPq) = |1) (S> |1)) the classical phase 
space distribution, given by the Husimi distribution, is 

-<-°)-S^nH^ • (30) 

k=i ^ ^ 



In the analytical calculation of the CSLMI we use the super-operator method 
to conform with classical Poisson brackets computations. We find 



modified 



/(t) = 8M(t)[2-8M(t)] (31) 

and 

r\t) = um - u{f)\ (32) 

where 

nit) = ^^^^^^ [5 + 3 cos (4At)] = 5f (t) = Sf{t) . (33) 

The quantum and the classical statistical linear entropies have similar qualitative behav- 
iors, but the amplitude of their oscillations are markedly different. They exhibit the same 
purification period tq = ^• 

Finally, we consider one of the subsystems in a coherent state and the other in the number 
state (I'i/'o) = ® |1)). In this case, for ai = 0, we obtain 



sin^(2At) .„ / , , XT 
I{t) = ^[7 + cos(4At)] 



, sin^(2At) , / , M 

P\t) = ^— ^ [15 + cos (4At)] 



(34) 



64 

These quantities are plotted in Fig. (jH). The qualitative agreement of the oscillations is 
remarkable, though the amplitudes are quite different. If the classical information is re- 
scaled so that its maximum at it/A coincides with the quantum plot, the two curves become 
very similar. This might seem a consequence of the normalization introduced in and not 
present in I (see Eq. (|2H) ). Unfortunately this is not so: the classical normalization is indeed 
necessary (at least for dimensional reasons) and choosing it in a case by case basis does 
not seem to bring any important physical information. Our interpretation of the differences 
in the amplitudes is that they reflect the non-classical character of the initial phase space 
distribution. 
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IV. CONCLUSION 



In this work we have defined the classical statistical linear mutual information, a tool 
that quantifies the non-separabihty of classical statistical distributions representing pure 
states of a bipartite Hamiltonian system. The comparison of the quantum and classical 
linear mutual information provides a measure of how much of the quantum loss of purity are 
due to intrinsic quantum effects and how much is related only to the probabilistic character 
of the initial distributions. We computed the classical and quantum mutual information 
for a system of two oscillators subjected to different types of coupling. We found that the 
two measures follow each other closely in the case of initially separable Gaussian states. 
For the case of linear coupling the classical and quantum mutual information are identical, 
revealing the classical nature of the system and the coherent evolution of the Gaussian wave- 
packets. For non-linear couplings the classical mutual information follows the quantum one 
for short times. This is follows from the fact that the short time quantum evolution can 
to.n,ulated . ... of L.ouviUe fo..a,.. Q Q. TM. p.o™jty . de..able fo. a 
measure of classical correlations in view of Ehrenfest's theorem 27[, and confirms that 
our definition is appropriate. For longer times the folding of the wave-packets certainly 
introduces self- interferences that have no classical counterpart. When quantum interferences 
become substantial ,24] the classical and quantum mutual information become significantly 
different. For open systems, where interferences are eliminated by the coupling with an 
external environment, we conjecture that the classical and quantum mutual information are 
going to coincide for much longer times. 

We have also investigated the role of non-Gaussian types of initial distributions in the 
classical separability. We have shown that also in this case the CSLMl Ici(t) is a meaningful 
quantity to measure the classical non-separability. Our example of linearly coupled har- 
monic oscillators (RWA) shows that the time evolution of the classical mutual information 
is qualitatively similar to that of its quantum counterpart. The amplitude of the classical os- 
cillations, however, is markedly different from the quantum ones, reflecting the non-classical 
nature of the initial state. 
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FIG. 1: QLMI I{t) and CSLMI Ici{t) as a function of time (arbitrary units) for the linear coupling 
and initial Gaussian distribution. The two curves are identical. 
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FIG. 2: Poincare section for the Nelson potential at = 0.05, gi = and > 0. The symbols 
stand for the centers of the coherent states: filled circle for chaotic region and filled square for 
regular region. 
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FIG. 3: QLMI (continuous line, h = 0.05) and CSLMI (dashed line) for the nonlinear coupling 

and initial Gaussian distributions centered on a chaotic orbit (a) and on a regular orbit (b). 
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FIG. 4: Quantum and classical mutual linear information for the RWA coupling for w = A = ^ = 1 
and initial state \iJjo) = \0)<Si\l). The continuous line shows the quantum result I{t) and the dashed 
line the classical P^t). 
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